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We check the eigenvalue spectrum of the 3>f +1 Hamiltonian against Poisson or Wigner behavior predicted 
from random matrix theory. We discuss random matrix theory as a tool to discriminate the validity of a model 
Hamiltonian compared to an analytically solvable Hamiltonian or experimental data. 



1. Motivation 

The fluctuation properties of the eigenvalues 
of a Hamiltonian give much insight into the dy- 
namics of the underlying system. In particular, 
the so-called nearest-neighbor spacing distribu- 
tion P(s), i.e., the distribution of spacings s be- 
tween adjacent eigenvalues plays an important 
role. According to the Bohigas-Giannoni-Schmit 
conjecture [1], quantum systems whose classical 
counterparts are chaotic have a P(s) given by 
random matrix theory (RMT) whereas systems 
whose classical counterparts are integrable obey 
a Poisson distribution, P(s) = e~ s . In this sense, 
the form of P(s) characterizes quantum chaos. 

Today we know a number of physical sys- 
tems where Wigner and Poisson behavior is ob- 
served [2]. Here we test a one-dimensional chain 
of N s coupled harmonic oscillators, with anhar- 
monic perturbation [3]. Its Euclidean action is 
given by 
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In the continuum formulation it corresponds to 
the scalar model, 
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Introducing a space-time lattice with lattice spac- 
ing a s and a t , this action becomes 

N 3 N t -1 

S = ^ ^ a t a s 

n=l k=Q 

1 / <$>(x n ,t k+1 ) - <f>(x n ,t k ) 

2 V a t 
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2 

+ ^ 2 ( Xrn t k ) + ^\x n ,t k ) . (3) 

The actions given by Eq.(3) and Eq.(l) can be 
identified by posing <f> — yfa^Q, Q — l/a s , f^o = 
m, and A/2 = g/4!. 

For the underlying real and symmetric matrix 
one expects a correspondence to the Gaussian or- 
thogonal ensemble (GOE). The RMT result for 
P(s) is quite complicated; it can be expressed 
in terms of so-called prolate spheroidal functions, 
see Ref. [4] where P(s) has also been tabulated. 
A very good approximation to P(s) is provided 
by 



P(s) 
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which is the Wigner surmise for the GOE. 

We computed the eigenvalues of the Hamilto- 
nian in Eq. (3) via the Monte Carlo Hamiltonian 
method using a stochastic basis [5] . This is a con- 
cept to calculate eigenvalues and eigenstates (in 
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Figure 1. Nearest-neighbor spacing distribution P(s) for a chain with N s = 9 harmonic oscillators from 
the analytical spectrum (left plot) and from the stochastic method (right plot). The Poisson and the 
Wigner distribution for the GOE are inserted. 
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Figure 2. P{s) for a chain with A^ s = 3, 5, 7, 9 anharmonic oscillators with A = 1 using a stochastic basis. 
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some energy window) of field theories or many- 
body systems. 

To construct P(s), one first has to "unfold" the 
spectra [1]. This procedure is a local rescaling of 
the energy scale so that the mean level spacing is 
equal to unity on the unfolded scale. Ensemble 
and spectral averages (the latter is possible be- 
cause of the spectral ergodicity property of RMT) 
are only meaningful after unfolding. 

2. Results 

Figure 1 shows the nearest-neighbor spacing 
distribution P(s) for the chain of harmonic oscil- 
lators with A = 0. This case reduces to the Klein- 
Gordon model which is classically integrable and 
thus a free theory. The quantized system of N s 
oscillators has a highly degenerate spectrum. We 
keep in mind that a single harmonic oscillator is 
an exceptional situation which would lead to a 
5-function at s = 1 after formal unfolding. The 
left plot in Fig. 1 depicts P(s) of the analytically 
known spectrum for N s = 9 oscillators. It might 
consist of a washed out peak from the harmonic 
oscillator spacing and of a peak around the origin 
from the Poisson distribution of the lifted degen- 
eracies. The right plot in Fig. 1 presents P(s) 
for the corresponding spectrum obtained from 
the stochastic basis (breaking parity and transla- 
tional invariance) with 1000 elements. It clearly 
exhibits Wigner behavior with the underlying 
symmetry of the GOE. It seems that the "holes" 
in the matrix produced by the method give ef- 
fectively a random matrix. Although the individ- 
ual exact eigenvalues are well reproduced [5] their 
fluctuations are different leading to a completely 
different dynamics. 

Figure 2 shows the nearest-neighbor spacing 
distribution P(s) for the oscillator chain when the 
anharmonic term is switched on, A = 1. We vary 
the length of the chain from N s = 3 to N s = 9 
and take 100 — 200 states for the stochastic ba- 
sis (with parity symmetry implemented) into ac- 
count. The plots are inconclusive and might show 
a trend to uncorrelated eigenvalues compared to 
A = which could be an effect from the deeper 
quartic potential. It will be interesting to see if 
this is enhanced with increased A. 



In summary, we have analyzed the spectrum 
of a finite chain of oscillators, with and without 
a quartic coupling. The infinitely long chain of 
harmonic oscillators is equivalent to the Klein- 
Gordon model being a free theory. Its spectrum is 
known analytically also for finite N s . We tested a 
stochastic method and obtained a spectrum cor- 
responding to a Wigner distribution in contrast 
to the theoretical expectation. Only preliminary 
results could be collected for the anharmonic os- 
cillators and conclusions have to be drawn in the 
future increasing the stochastic basis and the cou- 
pling space. 

This study is part of a general concept where 
random matrix theory is used to distinguish be- 
tween the quality of a model compared to a theory 
or experiment. An analysis in this respect is be- 
ing performed with a quark potential model and 
the experimental hadron spectrum [6]. 

Acknowledgments: This work was supported by 
NSERC Canada (H.K.), NNSF China (X.Q.L.j, 
and by FWF project P14435-PHY (R.P.). We 
thank G. Akemann for helpful discussions. 

REFERENCES 

1. O. Bohigas, M.J. Giannoni, and C. Schmit, 
Phys. Rev. Lett. 52 (1984) 1; O. Bohigas 
and M.J. Giannoni, Lcct. Notes Phys. 209 
(Springer, Heidelberg, 1984) 1. 

2. T. Guhr, A. Miiller-Groeling, and H.A. Wei- 
dcnmiillcr, Phys. Rcpt. 299 (1998) 189; E. 
Bittner, H. Markum, and R. Pullirsch, hep- 
lat/0110222. 

3. T.H. Scligman, J.J.M. Verbaarschot, and 
M.R. Zirnbauer, Phys. Rev. Lett. 53 (1984) 
215; J. Phys. A: Math. Gen. 18 (1985) 2751. 

4. M.L. Mehta, Random Matrices, 2nd ed. (Aca- 
demic Press, San Diego, 1991). 

5. C.Q. Huang, H. Kroger, X.Q. Luo, and 
K.J.M. Moriarty, Phys. Lett. A 299 (2002) 
483. 

6. H. Markum, W. Plcssas, R. Pullirsch, B. 
Sengl, and R.F. Wagenbrunn, in preparation. 



